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report was submitted 4 October 1966 by the AFWL Project Officer, Mr. Charles E. 
Needham (WLRTH). 
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ABSTRACT 


The results of a theoreticai calculation of the detonation 
of 1,000 pounds of TNT (loading density 1.608 gms/cc) are 


presented. The charge was detonated 15 feet above ground 


with an ambient pressure of 12.6 psi and an ambient tempera- 


ture of appproximately 100°F. The calculation started with 
the burning of the TNT and was carried to 115 milliseconds. 
The calculation clearly shows Mach stem formation, triple 


point path, and flow field. The theoretical calculation 


agrees well with experimental data obtained frcem a test 
conducted by AFWL at Holloman AFB. 
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SECTION I 


INTRODUCTION 


In recent years much effort has been expended in attempts to experimentally 
measure the parameters of the Mach stem and associated phenomena. Most recently 


this was done in Operation Distant Plain* (reference 1). 


In support of this operation the Air Force Weapons Laboratory was requested 
by the Defense Atomic Support Agency to make predictions of air blast parameters 
and associated two-dimensional errects to assist in pre-shot planning. These 
predictions are based oa theoretical calcuistions made by large hydrodynamic 
computer codes, which begin from first principles rather than from scaling laws 


or empirical models. 


After the theoretical calculations of Distant Plain phenomenology ‘ere com- 
pleted, Colonel Ralph Pennington learned of a similar experiment just completed 
at Holloman AFB, New Mexico. In this experiment a 1,000-pound spherical charge . 
of TNT was detonated and extremely fine field data were obtained. Therefore 
Colonel Peanington suggested a calculation be made of the Holloman shot for pur- 


poses of comparing theoretical calculations and experimental data. 


Presented in this report are the results of and comparison between the 


theoretical calcu'’ations and the experimenta) data of the Zollocan cxperiment. 
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® 
At this time, no experimertal data have been received from the Distant Plain 
series. 


ce eae wie’ 


eat PO aMMUON BE: ts haiti Nellneemmardnyrionie 


se erento eee Ten ree abe eenminaaaae eBaeM cet wee 


AFWL-TR-66-128 


SECTICN I1 


THE ONE-DIMENSIONAL CALCULATION 


The one-dimenaton=! calculation was made using 4 modified version of the 
AFWL SAP Code on the AFWL CDC-6600 computer. SAP is a one-dimensional, spherically 
symmetric, Lagrangian hydrodynamic computer code. It has been modified to include 


two asterials and a burn routine. 


For this calculation 800 zones were given equal delta R's (where R is the 
radius) of 0.16 cm. The first 254 zones were given a density of 1.608 gm/cm? and 
were flagged as TNT zones. These conditions resulted in a charge weight of 996 
pounds, the same weight as that used in the experiment. The remaining zones were 


given ambient atrospheric conditions. 


Although a hydrodynamically stable real atmosphere is a standard part of the 
SAP program, it was not possible to make use of this because of the unusual condi- 
tions found at Holloman AFB. The low pressure and hign temperature of the air 
at Holloman correspond--through the SAP equation of state for air--to a density 
which is far lower than that in a normal sea-level-based atmosphere for che same 
altitude (7000 feet). A stable atmosphere, which coincided with the 12.6 psi and 


100°F conditions on the day of the experiment, was developed. 


To calculate the burning of TNT, two input numbers were required: (1) the 
velocity of the detonation wave (a constant determined by density) in TNT, (2) 
the energy released per graew of TNT upon burning. The equatiosc cf state for TNT 


was als. required. 


The detonation velocity used for this run was 6.98 * 10° cm/sec, and the 
energy released was 4.26 « 10!= ergs/gz. The equation of state used was the LSZK 
formulation for TNT. (References 2 and 5.) For sa detailed description of the 


SAP code ind of the burr routine see reference 4. 


The detonation wave reached the surface of the TNT art 0.058 esec with a peak 
pressure of 1.75 ° 10° dynes/ca. Figure 1 shows the wave rr-file just before 
{it reaches the TNT surface, and figure 2 sin «s the profile after the front breaks 
the surface cf the TNT. In this figure the velocity is seen to increase in this 
free expansion phase. The pressure and density curves have the same shape that 


they did pricr to free expansion with the excepticn of the peak, which has been 
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eaten away by the rarefaction wave. The free expansion cools the INT and the TNT 
behaves like a cold piston compressing the low density air ahead of it. Figure 3 


shows the prefiles just before the rarefaction wave reaches the center. 


Fi,ure 4 shows the profile at the beginning of the two-dimensional calcula- 
tion. Notice that all velocities are positive and that a sharp velocity gradient 
has developed just inside the TNT-Air interface. 


The one-dimensional calculation was continucd beyond the beginning of the 
two-dimensiona] calculativn to compare the two runs and to show more clearly the 
two~dimensional effects. Figure 5 shows the profile as a ehock moves inward from 
the TNT surface to the center. This shock is reflected from the center and forms 
a second outgoing shock. Many more small shocks are formed and follow in a simi- 
lar manner. Following the reflected shocks in the one-dimensional case is 
interesting but probably of little value since two-dimensional effects are much 


more important. 
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Figure 3. Profiles of Hydrodynamic Variables as a Function of Distance 
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i SECTION III 


y THE TWO-DIMENSIONAL CALCULATIC 4 


A The two-dimensional calculation was made using the APWL SHELL-OIL code on a 

- CDC-6600 computer. SHELL-OIL is a two dinensional, pure Eulerian, axially sya- 

metric, hydrodynamic code. SHELL is a one-material code; that is, ail material : 
in the problem is assumed to have the same equation of state. For this calculation 
the Poan-Nickel equation of state for air (reference 5), an empirical fit to 


Hilsenrath's data, was used. (For a detailed discussion of SHELL see reference 4.) 


Some question arises as tc the validicy of treating TNT as air. At the tine 
when the two-dimensional calculation was started (0.95 msec), the TNT had expanded 
to such an extent that it had a density comparable to that of the ambient air. 
This fact, combined with the relatively low temperature found thrcughout the TNT 
at this time. insures the validity of this assumption. To confirm this, a short 
one-dimensional calculation (SAP) was made using the air equation of state for all 
material. Differences between the two one-dimensional runs were on the order of 


l percent. 


A minor modification was made in the SHELL program to allow monitoring at 
fixed points in the grid. Nineteen such points (test stations) were used in this 
calculation. The first eleven test stations were chosen to coincide with the 
experimental instrumentation stations. The twelfth through nineteenth test 


Stations were placed at burst height at increasing radii. A coeplete list of 


station positions is given in tabie I. 


The plets from these stations show the initial shack, the shock reflected 
from the ground, the second shock from the center, and the passage of the triple 
point past burst height. Figure 6 shows the triple point path relative to the 


test stations. 


Included ir this calculation were 801 trace particles. These particles fol- 
low the fluid mot{on but do not influence the hydrodynamics in any way. Initially, 
these particles were placed on the TNT-Air interface at equal intervals. The sove- 
ment of these particles represents the motion of the interface as a function of 
time and appears as a heavy line or the contour and velocity vector plots in 


? 


Appendix I. 
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Table I 


TEST STATION POSITIONS AL 


Ground Range Height | 
Station Meters) | Reet. be -sGisteray | Feet 
i (10.67) 35.0 (1) 3.28 
2 (12.8) 42.0 (1) 3.28 ae 
3 (15.5) 51.0 (1) 3.28 fo 
4 (17.98) 59.0 (1) 3.28 | 
5 (19.5) 64.0 (1) 3.28 
6 (21.03) 69.0 (1) 3.28 
7 (23.47) 77.0 (1) 3.28 
8 (26.52) 87.0 (1) 3.28 
9 (29.26) 96.0 (1) 3.28 
10 (31.39) 103.0 (1) 3.28 
ll (35.35) 116.0 (1) 3.28 
12 (5.0) 16.40 (4.57) 15.0 
13 (6.0) 19.69 (4.57) 15.0 ; 
14 (8.0) 26.25 (4.57) 15.0 
15 (10.9) 32.80 (4.57) 15.0 : 
lo (12.0) 34.37 (4.57) 15.0 
17 (15.0) 49.21 (4.57) 15.0 
18 (20.0) 65.62 (4.57) 15.0 
19 (25.0) 82.02 (4.57) 15.0 
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Initially the SHELL mesh was a rectangle 142 zones in the radial (r) 
direction by 156 zones in the vertical (z) direction, in which each zone was 
6 co square. The entire rectangle has cylindrical symmetry about the verticle 
line r= 0. A sphere of radius 3.5515 m centered at r = 0, 2 = 4.572 m was 
generated. Into this sphere were placed the hydrodynamic variables as a function 
of radius, as calculated by SAP at 0.95 msec. Outside of this sphere the same 
ambient atmospheric conditions were generated as had been used previous}y in SAP. 
The bottom boundary, 2 = 0, was made reflective, and the calculations was begun. 


As the calculation progressed and the shock expanded, an automatic rezone 
enlarged the grid by approximately 6-1/2 percent in linear extent each time the 
shock reached the edge (top or side) of the grid. 


The saall periodic jumps in the overpressure versus time curves and the 
associated plots in Appendixes I and II were due to this rezoning. The jumps 
occurred when a station, fixed in space, changed its position relative to the 
mesh by an amount sufficient for the station to occupy a different zone. This 
condition, in effect, changed the position of the station by one zone dimension. 
The jumps could have been smoothed very easily but would have only served aesthetic 


purposes. 


In the discussion of the one-dimensional calculation, a shock was shown 
converging on the center of burst and forming a second outgoing shock. The radius 
of the outer shock was approximately 30 feet st this time. This means that the 
outer shock reached the ground and reflected to the center of burst before the 
converging shock had tide to reach the center. Consequently, the convergence 
occurred at a point higher than the original center and alec earlier than in the 
one-dimensional case. However, the second shock having been reflected from this 
new center, reached the ground at che same time in both calculations. It is this 
second shock that appears ss the emaller peak on the overpressure versus time 
curves for stations 1-11. (Appendix II.) 
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SECTION IV 


COMPARISON OF EXPESIMENT AND CALCULATION 


The experiment was conducted by AFWL ia conjunction with the Blast Environ- 


went Development portion of the Rocket Sled Blast Simulation Program at Holloman 
AFB, New Mexico. 


A 996-pound sphere of INT was detonated 15 feet above ground levei. The data 
collected include the high-speed photographs and pressure-time histories at the 
first 11 stations listed in table I. All experimental instrumentation was below 
the triple point and thus in the Mach stem. 


The analog data tapes were digitized, and plots were made at the Kirtland AFB 
Date Reduction Center. Arrival times and positive phase durations were read 
directly trom the digitized tapes, but no direct measurements of overpressure 
impulse were mac>. The impulses given on the data graphs were computed by numeri- 
cal integration of the overpressure cime curves. Because of instrument noise, 
these impulses had oscillations of about 20 percent and have been included only 
for completeness. 


Figure 7 shows the peak overpressure, both experimental and calculated, as 

a function of ground range. The overpressure impulse is that calculated by SHELL. 
The maximum difference between experimental and calculated overpressure is approxi- 
mately 18 percent. The experimental curve ceciliates about the calculated curve i 
and the resulting differences have been interpreted as exper:mental error. Teble 

II contains the exact values for overpressure and overpressure impulse data. Figure i 
8 shows ihe arrival time and positive phase durat/on of the mgin shock, both experi- 
wental and calculated, as a function of yround range. Again, the experisental ; 
data oscillate about the calculated curve with a saximue deviation of 7 percent ir 


arrival tiwe and of 12 percent in positive phase duration. The difference in 


arrivel time is within experimental error. The larger percentage error in poa‘tive 
phase duration hae beer attributed to che compounded error jn ar:ival time and to 
the electronic ruise in the overpressure seasurements. The roise oresent  sused 


the ow: cspressure to appear uegative in some instances before and others after the 
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Table It 


OVERPRESSURE AND OVERPRESSURE IMPULSE DATA 


Distance : 
Station _ Pasoce | tare ——_J op oP: 


1 35 98.2 153.4 95.8 167.0 
2 42 61.5 138.8 63.8 151.0 
3 51 40.7 115.9 41.9 128.4 
4 59 32.5 117.2 29.2 111.9 
5 64 21.3 91.4 24.2 104.8 
6 69 18.2 90.4 29.8 96.0 
- 7 77 19.8 105.5 16.9 89.8 
8 87 13.1 $7.9 13.2 80.2 
9 96 9.5 69.9 12.0 74.5 


~ 
i=) 


103 8.1 69.0 9.6 67.7 
116 7.4 38.0 7.7 62.8 
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The depression at the back of the shock on the experimental curves was caused 
by instrument design aad should be ignored. This depression, however, contributed, 
dy lowering the impulses, to the previously described error, in the experimental 
overpressure impulse calculation. Figure 9 shows dynamic pressure and dynamic 


pressure impulse veraus ground range. 


The greatest discrepancy between the SHELL calculation and the experiment 
is due to the second shock arrival time. The calculation indicates a second 
shock arrival time about 3 milliseconds greater than thet measured experimentally 
at the first station, and this difference in time increases with distance. (Com- 
pare the experimental and theoretical graphs in Appendix II.) At this time, no 
satisfactory explanation has been made for the difference, but several possibili~ 


Ra ties are being investigated. 


The minimum overpressure before second shock arrival and the peak over- 
pressure of the second shock agree within experimental error. This makes the 
ee discrepancy in second shock arrival time even more diffiiilt to understand and 
oe explain. 
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SECTION V 


CONCLUSIONS 


z. The experimental data included in this report are believed by the experi- 
menters to be some of the best data of this type available. If this is the case, 
one can conclude that the parameters measured in this type of experiment are more k 
easily and accurately determined by SHELL than by experimental measurement. 
Further, parameters not measured in this experiment, such as dynamic pressure, 
dynamic pressure impulse, and velocity, are also more easily and precisely deter- 
mined by SHELL. Even though the discrepancy in the second shock arrival time is 
of secondary significance, it should not be ignored. This discrepancy does not 


\ affect the principal parameters of the problem. 


Included in Appendix I are photographs of the experiment and calculated con- 
tours both occurring at approximately the same times. The initial shock reflec- 
tion, Mach stem formation, and second shock reflection are much more easily 


distinguished on the contours than on the films. 


In view of the comparisons between the SHELL wethod and experimental measure- 
ment presented in this report, much time, effort, and expense could be saved with 
pre-experiment calculations. Instrument calibration and position could be deter- 
mined to achieve maximum utilization. By using the results of a calculation such 
as this, much better field data could be collectcd, and from these field data the 


codes could be improved to make even more accurate calculations. 
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Appeudix I 


CONTOURS, VECTORS, AND PHOTOS 


The first three figures contai:.ed in this appendix are photographs of the 
detonation at 1, 2, and 3 milliseconds. The edge of the fireball in the photo- 
graphs corresponds to the edge of the TNT detonation products. These in turn 
correspond to the particles on tne interface seen in the computer plots, but the 
air shock is not visible in these photographs. Comparisons between photographs 
and plots are possible but diifzicult due to the qua*\ity of the film. A more use~ 
ful comparison can be fou: 1 in Appendix II. 


The plots included in Appendix I indicate the formation and growth of the 


Mach stem with th2 associated velocity, pressure, and density distributions. 
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Appendix ITI 


STATION 1 TO 11 PLOTS 


This appendix includes the experimental overpressure vs. time curves at the 
eleven test stations and the computed curves taken at the same positions as at 
the experimental stations. Computed curves include overpressure, two components 


of dynamic pressure, the corresponding impulses, and the two components of velocity. 


This appendix presents the most consistent comparison of experiment and cal- 


culation. A summary of this comparison is found in figures 6 and 7 of the text. 
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Appendix III 


STATION 12 TO 19 PLOTS 


The stations uced to obtain these plots were located at burst height with 
increasing radii. 
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